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ANALYSIS OF A COUPLED SPIN DRIFT-DIFFUSION 
MAXWELL-LANDAU-LIFSHITZ SYSTEM 


NICOLA ZAMPONI AND ANSGAR JUNGEL 

Abstract. The existence of global weak solutions to a coupled spin drift-diffusion and 
Maxwell-Landau-Lifshitz system is proved. The equations are considered in a two-dimen¬ 
sional magnetic layer structure and are supplemented with Dirichlet-Neumann boundary 
conditions. The spin drift-diffusion model for the charge density and spin density vector 
is the diffusion limit of a spinorial Boltzmann equation for a vanishing spin polarization 
constant. The Maxwell-Landau-Lifshitz system consists of the time-dependent Maxwell 
equations for the electric and magnetic fields and of the Landau-Lifshitz-Gilbert equation 
for the local magnetization, involving the interaction between magnetization and spin den¬ 
sity vector. The existence proof is based on a regularization procedure, L^-type estimates, 
and Moser-type iterations which yield the boundedness of the charge and spin densities. 
Furthermore, the free energy is shown to be nonincreasing in time if the magnetization- 
spin interaction constant in the Ladau-Lifshitz equation is sufficiently small. 


1. Introduction 

Magnetic devices, such as magnetic sensors and hard disk read heads, typically consist of 
ferromagnetic/nonmagnetic layer structures. A model for magnetic multi-layers was hrst 
introduced by Slonczewski [33]. This model is well suited for Magnetoresistive Random 
Access Memory (MRAM) devices but it is less appropriate for current-driven domain wall- 
motion. A more general approach is to introduce the spin accumulation coupled to the 
magnetization dynamics. The evolution of the magnetization is modeled by the Landau- 
Lifshitz (-Gilbert) equation [3H]. When electrodynamic effects cannot be neglected (like in 
high-frequency regimes), this description needs to be coupled to the Maxwell equations. 
In this paper, we analyze for the hrst time a coupled spin drift-diffusion Maxwell-Landau- 
Lifshitz system in two space dimensions with physically motivated boundary conditions. 

Let us describe our model in more detail. We consider a three-layer semiconductor 
structure C consisting of two ferromagnetic regions Wi, LjJ2 C fl, separated by a 
nonmagnetic interlayer f2\ci;, where cn = cui U Ci ;2 is the union of magnetic layers [1]. 
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Landau-Lifshitz-Gilbert equation. The dynamics of the magnetization m = (mi, 
m 2 , m 3 ) is governed by the Landau-Lifshitz-Gilbert (LLG) equation 

(1) dtm = m X (Am -|- H -|- /3s) — am x (m x (Am -(- H + /3s)) in cu, t > 0, 

where the effective held Hefj = Am -|- H consists of the sum of the exchange held contri¬ 
bution Am and the magnetic held H, and a > 0 denotes the Gilbert damping constant. 
The additional term /3s models the interaction between the magnetization m and spin 
accumulation s with strength /3 > 0 [9l[36]- We choose the initial and boundary conditions 

(2) m(0) = m° in a;, Vm ■ 1 / = 0 on u, t > 0, 

where jz is the outward unit normal on du, we write m( 0 ) = m(-, 0 ), and the notation 
Vm ■ u = 0 means that Vm* ■ u = 0 for 3 = 1,2, 3. The Neumann conditions were also 
used in, e.g., mm- We set m = 0 in G\a;. 

The existence and non-uniqueness of weak solutions to the LLG equation goes back to 
[31 [31]. The local existence of a unique strong solution was proven in [5]. In two space 
dimensions and for sufficiently small initial data, the strong solution is, in fact, global in 
time |1|. For general initial data, the two-dimensional solution may develop hnitely many 
point singularities after hnite time; see [20] for a discussion. The existence of weak solutions 
in three space dimensions with physically motivated boundary conditions was shown in [ 1 ], 
based on a hnite-element approximation. For a complete review on analytical results, we 
refer to [TOl [26] . 

Maxwell equations. The Maxwell equations are given by the time-dependent Ampere 
and Faraday laws for the electric and magnetic helds E = {Ei, E 2 , E 3 ) and H = 

H 3 ), respectively, 

(3) — curl H = Jg, c3iH-|-curl E = — c3fm in G, t > 0, 
and by the Gauss laws 

(4) divE = p - G(a;), div(H-h m) = 0 in G, t > 0. 

Here, Jg is the electron current density, p the electron charge density, and C(x) the doping 
concentration characterizing the device under consideration. We assume that the boundary 
dfl splits into two parts: the Ohmic contacts F £, and the union Fat of the insulating parts, 
with dQ = F/) U Fat. Then the initial and boundary conditions of E and H are given by 

(5) E(0) = E°, H(0) = in G, 

( 6 ) E X 1 / = 0 on Fa), f > 0, H x i/ = 0 on Fat, t > 0, 

(7) E ■ 1 /= 0 on Fat, t > 0. 

The existence analysis (for given and smooth Jg and dtui) may be based on Kato’s theory 
of quasilinear evolution equations of hyperbolic type [2H] or on semigroup theory [22]; also 
see Section [ST] 

Goupled Maxwell and LLG equations were intensively studied in the literature. For 
instance, the Maxwell-Landau-Lifshitz system in three space dimensions with periodic 
boundary conditions was investigated in [19]. Garbou and Fabrie [ 6 ] proved the existence 
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of weak solutions to the LLG equation, coupled to Maxwell’s equations, in the whole space 
The existence of spatially periodic strong solutions in three dimensions and their local 
uniqueness were proved in [9]. The solutions are only partially regular (i.e. smooth except 
on a low-dimensional set) because of possible vortices or phase transitions. We refer to [I3] 
for the two-dimensional case and to [ElIIl] for three space dimensions. 

Spin drift-diffusion system. We consider the spin drift-diffusion equations for the 
charge density p and the spin density vector s = (si, S2, S3) 

(8) dtp-div Je = 0, Je = T)(Vp - pE), 

(9) dtS — div Js -I- 7m X s = —s/r, = T)(Vs — s ® E), 

where D > 0 is the diffusivity constant, J<j the spin current density vector, 7 > 0 is the 
strength of the effective magnetic held, and r > 0 denotes the spin-hip relaxation time. 
The term 7m x s causes the spin density vector to rotate around the magnetization, while 
the spin-hip relaxation term leads, in the absence of other forces, to exponential decay to 
the equilibrium spin density vector Seq = 0. We assume that the densities p and s are 
prescribed on T^i (Ohmic contacts), while there are no-hux boundary conditions on T^r 
(insulating boundary). This results in the initial and boundary conditions 


(10) 

1>(0) 

0 

0 

in fl. 


(11) 


p = Pd, s = 0 

on Td, 

f > 0, 

(12) 

T- 

u = 0, Js ■ = 0 

on Tat, 

t > 0. 


The spin current density is a 3 x 3 matrix with rows Js,i = Vs* — SjE for i = 1,2,3. 
Accordingly, ■ 1/ is a vector in consisting of the elements Js,i ■ = 0, i = 1, 2, 3. 

Spin-polarized drift-dihusion models were analyzed only recently in the literature. Glitz- 
ky [TS] proves the existence and uniqueness of weak solutions to a two-dimensional transient 
drift-dihusion system for spin-up and spin-down densities. The stationary problem was 
solved in three space dimensions in [TB]. These models were derived from the spinor 
Boltzmann equation in the dihusion limit with strong spin-orbit coupling in [15] . 

More detailed information can be obtained by introducing the spin density. Spin-vector 
drift-dihusion equations can be derived from the spinor Boltzmann equation by assuming a 
moderate spin-orbit coupling [15]. Projecting the spin-vector density in the direction of the 
magnetization, we recover the two-component drift-dihusion system as a special case. In 
[TB] . the scattering rates are supposed to be scalar quantities. Assuming that the scattering 
rates are positive dehnite Hermitian matrices, a more general matrix drift-dihusion model 
was derived in [30]. The global existence of weak solutions to this model was shown in 
[25] . An energy-dissipative hnite-volume discretization was presented in [7]. 

Equations (j8D-([9|) result from the cross-dihusion model in [30] by choosing a vanishing 
spin polarization constant. By this choice, the dihusion matrix becomes diagonal which 
makes our analysis possible. For a more general spin drift-dihusion LLG model, but without 
coupling to Maxwell’s equations and with saturating drift velocity, we refer to [35] . 
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In the physical literature, also other equations for the spin density vector have been 
suggested. In [2S1 Formula (8)], the spin density is dehned as the difference of the spin- 
up and spin-down densities. Thus, the underlying equation is a two-component model 
which is a special case of the general model. Starting from kinetic equations for the 
charge and spin components of the Wigner-transformed density matrix, Lueffe et al. m 
Formula (54)] derived a spin diffusion equation for weak spin-orbit interaction or strong 
scattering. Another derivation employs a SU (2) gauge held theoretical description of the 
spin-orbit coupling and the Heisenberg held operators for the dehnition of the spin density 
[SU Formulas (l)-(4)]. The resulting equation is similar to ([9]) but the spin current density 
also depends on the charge current. Finally, assuming that the dihusivity in the drift- 
dihusion equation for the density matrix is proportional to the magnetization vector, the 
authors in [SB] obtain ([9]) with a spin current density whose drift term equals m E 

instead of s (8) E as in ([9]) . The former drift term can be derived from the Wigner equation 

in the dihusion limit by approximating the Wigner function appropriately [H Formula 
(23)]. We stress the fact that the model (l8])-([9]) is derived from the spinor Boltzmann 
equation without heuristic arguments. 

Main results. We show that there exists a global-in-time weak solution to the coupled 
spin drift-diffusion Maxwell-LLG system. Our assumptions are as follows: 

(13) cu C O C are bounded domains with smooth boundaries, 

(14) a, /3, 7 , H, r > 0, C e T“(0), 

(15) Pd, s°, E°, G Fr^(O), m° G |m°| = 1 in oo, 

(16) divE° = — 0(x), div(H° -|- m°) = 0 in O, E° ■ = 0 on Fjv. 

We also suppose that dfl = F^iUF n, FdOFat = 0, and Ftv is open and has positive measure 
in dQ. To simplify the notation, we write sometimes u G H^{Q) instead of u G for 

vector-valued functions u. We denote by the space of all functions in H^{Q) with 

zero trace on To and by H\){VL)' its dual space. 

Let us discuss assumptions (lT3|) - (fT6D . The restriction to two space dimensions is (only) 
needed in the uniqueness proof for the regularized LLG equation fl25D . This property is 
required to obtain a well-dehned hxed-point operator. In (ITBD . we suppose that equations 
(jl]) and ([7]) hold initially. These properties allow us to conclude the validity of (jl]) and ([7]) 
from ([3]) and ([6]) (see e.g. m p. 435f.]). 

The first main result is the following theorem. 

Theorem 1 (Existence of global weak solutions). Let assumptions flTBD - flTBD hold. Then 
there exists a weak solution to CD-da satisfying 

p > 0 in Q, m = 0 in Q\u, t > 0 , 

p, s G LL(0, oo ; H\n)) G L;^,(0, oo ; L-(G)), d^p, d^s G 11^(0, oo; i7i(G)'), 

E, H G G°([0,cx));L2(G)), 

m G T“^(0, oo; dtm G oo; L‘^{u)), |m| = 1 m w. 
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The L°° bounds on p and s can be shown to be uniform in time; see Remark [HI 
The proof of this theorem is based on a combination of semigroup techniques for the 
Maxwell equations (I3j), a Galerkin approximation for the LLG equation ([I]), and esti¬ 
mates for the spin drift-diffusion model (E])-©. Note that it is sufficient to solve ([3]) with 
dS])-© as dl]) with dZ]) are consequences of the former equations. Since Je and dtva. are not 
regular a priori, we approximate these terms by regularizing Vp, dtm and truncating p, s 
in the drift terms in ©-©, respectively. This regularization is similar to that employed by 
Jochmann [22] for a coupled Maxwell drift-diffusion system (without spin). The challenge 
in the proof is to remove the regularization and truncation. For the de-regularization limit, 
we derive uniform estimates for the variables by showing that the functional 

(17) = ^ f [{p - pdY + \s\‘^ + \E\‘^ + \'}i\‘^)dx + I- f \Vm\‘^dx 

^ Jn ^ Jui 

satishes the inequality 

S'(t)-l-ci [ / (|Vpp-h |Vsp)da:-1-C 2 [ f \dtm\'^dx < Cs{T), tG(0,T), 

Jo Jn Jo Ju) 

where Ci, C 2 , c^{T) > 0 are some constants which are independent of the solution. Further 
details on the proof are given in Section [2l In order to remove the truncation, we derive 
L°° estimates for m, p, and s by using a Moser-type iteration procedure. 

The functional S{t) is not the energy of the system. The (relative) free energy consists of 
the von-Neumann energy for the spin system, the electromagnetic energy, and the exchange 
energy of the magnetization: 

(18) E{t) = ^f (p+(logp+- 1)-Fp_(logp_ - 1)-21ogpD(p-pD))dx 

+ ^ [ {\E-log pd\‘^ + \ii\^)dx+ ^ f \Vin\^dx, 

^ Jn ^ Ju) 

where p± = p ± |s| (see Section 0]). This formulation implicitly assumes that p > |s|. Our 
second main result is the proof that E[t) is nonincreasing in time under the conditions that 
the interaction parameter /? > 0 is sufficiently small and the solution is smooth and satishes 
p > |s|. This shows that the coupled system dissipates the free energy. The constraint 
on the parameter fj may come from the fact that the term fJs is introduced in the LLG 
equation only heuristically, but we leave further investigations to future research. 

The paper is organized as follows. The strategy of the existence proof is explained in 
Section [2] and the full proof is given in Section |3l We conclude in Section |4] with the 
monotonicity proof for the free energy E{t). 

2. Strategy of the proof of Theorem 0] 

In order to prove Theorem [H we hrst consider a truncated and regularized problem. For 
this, let T > 0, e > 0, M > 0 and set [z\m '■= min{M, max{0, z}} for 2 ; G M. We wish to 
prove the existence of weak solutions to 

(19) 9tp-div(T)(Vp-[ p]mE)) = 0, 
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(20) dts - div \^D 

(21) ^tE-curlH 

(22) ^iH + curlE 

(23) dtin — eAm = 


Vs 


|s|]M-y (2) E 


+ 7 m X 


s r 


- |p]mE), 

—dtRl{m) in 11 X (0, T),, 

m X (Am + H + /5s) — am x (m x (Am + H + /5s)) in w x (0, T) 


with the initial and boundary conditions ([2]), ([5])-([6]), and (fT0|) - ffT^ . In the Maxwell equa¬ 
tions (l2T]) - (l2^ . and R^ are two families of linear regularization operators acting on 
functions of x and t, respectively, satisfying for all u G L‘^{VL) and v G L^(0,T), 


(24) 

PeHIlcqn) < kelHlmn), 



(25) 


lim 

£—^0 

\\Rl{u) - ■■ 

(26) 

1 -^e(^)llci([0,T]) A fce||'?^||L2(0,T), 



(27) 

1 11^1(0,T) < ko\\v\\H^{0,T), 

lim 1 

£—^0 

1 1.2(0,T) 


where > 0 depends on e but fco > 0 is independent of e. The space-regularization 
operator R^ was introduced in [22l p. 665f], where also their existence and the above 
properties were proved. The time-regularization operator R\ can be dehned in a similar 
way. 

In the following, we abbreviate A = L2(0, T; L2(11))4 and Y = C^{[0,T]-, 

The first step of the proof of Theorem [T] is the application of the Leray-Schauder fixed- 
point theorem to the map 

F : A X y X [0,1] —)■ A X y, (p, s; m; a) i-A (p*, s*; m*), 

which is defined as follows (details will be given in the following subsections). Let (p, s; m; 
cr) G A X y X [0,1] be given. 

I. Solve the regularized Maxwell equations 


(28) dtE - curl H = aD(VR^(p) - Ip]mE), 

(29) dtH + curlE = —adtRl(m) in hi x (0,T), 
with initial and boundary conditions 

(30) E(0) = (tE°, H(0) = (tH° in Q, 

(31) Exi/ = 0 onT^i, f>0, Hxi/ = 0 onTjv, t>0. 


and obtain (E, H) G C°([0,T]; L^(Q)). 

II. Solve the regularized (nonlinear) LLG equations 

(32) dtm* — eAm* = m* x (Am* -|- aH -|- a/3s) 

— am* X (m* x (Am* -|- crH -|- cr/5s)) in a; x (0, T) 


with initial and boundary conditions 

m*(0) = crm° in cj, Vm* ■ 1 / = 0 on cj, t > 0, 


( 33 ) 
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and obtain m* G L°°(0, T; T ; satisfying |m*| < 1 in x (0, T). 

III. Solve the linearized spin drift-diffusion equations 


(34) 

(35) 


dtp* - dw{D{\/p* - a[p*]ME)) =0 in f] x (0, T), 


dts* - div { D { Vs* - o-[|s*nMT—7 (8 E 


cjym X [|s Wmttj = -> 


r 


with the initial and boundary conditions 


(Vp* 


cr[p*]ME) - v = D 


p*(0) = (jp°, s*(0) = crs° 

p* = crpD, s* = 0 

Vs* - (T[|s*|]Mj^ (8) E j ■ v> = 0 


in 


on Td, 

t > 0, 

on Td, 

t > 0, 


and obtain (p*, s*) e L‘^{0, T; H\n)) n T; ifi(fi)') C X. 

The regularization in (|2^ ensures that the solution is unique, which is necessary for 
the dehnition of the hxed-point operator. Observe that F(p, s; m; 0) = (0, 0; 0) since 
the solutions to the homogeneous subproblems (a = 0) are trivial. Standard arguments 
show that the operator F is continuous. By Aubin’s lemma |32], bounded sequences in 
L^(0, T; hf^(O)) n are relatively compact in L^(0, T; L^(f2)) and bounded 

sequences in L°°(0, T; niir^(0, T; L^(0)) are relatively compact in O°([0, T]; L^(f2)). 

Consequently, F is compact. It remains to prove uniform estimates for all hxed points of 
F{-,a). They will be derived from estimates for the functional S{t) dehned in flTT)) : see 
Section 13.31 Then the Leray-Schauder hxed-point theorem implies the existence of a hxed 
point of F(-, 1), i.e. of a solution to flT^ - fl^ with the corresponding initial and boundary 
conditions. 

The estimates from S{t) turn out to be independent of e which allows us in the second 
step of the proof to perform the limit e —)■ 0. The proof that we can remove the truncation 
in flT^ - fl2n]) is more delicate. We prove in the third step L°° bounds for p and s by employing 
a Moser-type iteration technique. The idea is to derive estimates of p and s, which are 
independent of p, and then to pass to the limit p ^ oo. By a rehned Moser-Alikakos 
iteration technique, it is even possible to show that the L°° bounds are uniform in time; 
see Remark [HI 

The proof of |m| = 1 is slightly diherent. First, we show for the regularized LLG 
equation ([23]), by using a Moser-type iteration, that < 1. After the limit e 0, 

we can take the inner product of the limit equation ([ID and m to deduce immediately that 
|m| = 1 in a;, t > 0. 


3. Proof of Theorem □ 

3.1. Analysis of the regularized Maxwell equations. We show that the regularized 
Maxwell equations f[?ID - f[22D are uniquely solvable. 
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Lemma 2 (Existence of the regularized Maxwell equations). Let (p, m) G L‘^{0,T] L‘^{Q)) 
xC°([0, T]; L^(r2))^ and a G [0,1]. Then there exists a unique mild solution (E, H) G 
C°([0,T; to (EHD-iEI]). 

Note that equations (jl]) and d?]) are yet not proved. They will be shown in Section 13.41 
to hold for the de-regularized system. 


Proof. The proof is based on semigroup theory and the Banach hxed-point theorem. In 
principle, a fixed-point argument is not necessary since the Maxwell equations are linear. 
However, we would need to deal with a non-autonomous operator because of the presence 
of the term [p(x,t)]ME on the right-hand side of (|28|) . Therefore, we prefer the simple 
hxed-point argument. Following [22], we introduce the spaces 

IT = {u G : curlu G 

We = |u G IT ; J u — u • curl (j))dx = 0 for 0 G Ar)^|, 

Wh = |v G it : y (^ ■ curlv — v ■ curlu)(ix = 0 for u G 1Te|. 

The space 1 Tb consists of all functions u satisfying uxi/ = 0onrBiina generalized sense, 
and the space Wh consists of all functions v such that vxiz = 0 onrjvina weak sense. 
It is shown in Theorem 1, Chapter IX, § 3 of [n] that the operator 

A : We X Wh x (u, v) i-A (— curl v, curlu). 


is skew self-adjoint, i.e. A* = —A. Thus, —iA is self-adjoint and by the Theorem of Stone, 
—A generates a unitary Cq group (e-*^)ieR in L2(fi)3 x 

The regularized Maxwell equations fl2SD - fl2^ can be reformulated as 

dt{E,H) + A{E,H) = a{D{VR^,{p) - [p]mE),- dtRl{ni)), t > 0. 


The right-hand side is a function in L^(0, T; L^(r2))®. Thus, by Duhamel’s formula. 


(E,H)(f) = e-'"^(E°,H°) + a / e 




(D(VWM - MME),-S,fl‘(m))(s)*. 


We infer that the solutions to are the hxed points of the operator G : C°([0,T]; 

L^{n)f C%[0,T];L‘^{n)f, dehned by 

G(E,H) = e-'^(E0,H0) + a T e-(*-*)^(D(Vi?,"(p) - [p]mE),- dtRl{in)){s)ds. 

Jo 

Since (e“*^)igR is a unitary group and 0 < [ p]m < M in n X (0,T), we obtain for (E, H), 
(E',H')GCO([0,T];L2(fi))6, 

rt 


||G(E,H) -G(E',H') 


1l2(0)6 — 


e-(t-M(/^[p]^(E-E'),0)(s)ds 


L 2 ( 0)6 


< DM / ||E - E'||B2(f^)3(is < T)MT||E - E'|lco([o,r];L2(Q))3. 
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Thus, choosing T > 0 sufficiently small, G becomes a contraction, and there exists a unique 
local-in-time mild solution (E, H) to fl25]) - fl2^ . The global solvability is a consequence of 
the energy estimate (see e.g. [211 Prop. 2.4]): 

~ [ (|E|2 + \u\^)dx = a [ (DiVRIip) - [p]mE) ■ E - dtRl{in) ■ B.)dx 
^ Jn Jn 

<a [ {DVR^,{p)-E-dtRlim) ■U)dx 

Jn 

<\l (|E|2+|Hp)* + c f (\S7R-(pp+\d,Rl{vn)f)dx, 

^ Jn Jn 

where here and in the following, c > 0 denotes a generic constant independent of e and 
M if not stated otherwise. By Gronwall’s lemma and the properties (l24j) and fl2^ of the 
regularization operators. 


(|E(f)p -F \'H.{t)\‘^)dx <c + c{e) 



{p^ + \m\^)dxds, f > 0. 


0 Jn 


This estimate allows us to continue the local solution for all time t > 0. 


□ 


3 . 2 . Analysis of the regularized LLG equation. We show that fl 32 l) possesses a unique 
strong solution. 

Lemma 3 (Existence of the regularized LLG equation). Let Hg. := (t(H -F jJs) G R^iuj x 
(0,T)). Then there exists a unique strong solution m to fl52]) - fl55]) satisfying |m| < 1 in 
cj X (0,T), 

m e L°°(0, T; H\uj)) G ^^(0, T; H^{uj)), dtui G L^(0, T; L^(u;)), 
and the estimate 

||m||L°°(0,T;J?l(a;)) + + \\dt™-\\L^{0,T-,L^{u))) < C, 

where c > 0 zs independent of e. 

Proof. The proof is based on the Galerkin method, standard Lf estimates, and a Moser- 
type iteration to prove the L°° bound for m. 

Step 1: Existence of solutions to fl32l) - fl3^ . Let e* G fl L°°{u) (z G M) be the 

eigenfunctions of —A in u with homogeneous Neumann boundary conditions and with 
associated eigenvalues A* > 0. Let m^^^ ix,t) = ej{x) be the approximated 

solution to that is 

(36) f — m^^^ x (Am^^^ -F Hg.) 

J (jJ 

-F am^^^ X (m*^^^ x (Am*^^^ -F iia-)))eidx = 0, i = 1,..., N. 

This is a system of ordinary differential equations in the unknowns tv, which 

has a unique solution m*^^) in a suitable time interval (0,T*) with T* < T. It remains 
to hnd A-independent estimates for m*^^) in order to conclude global solvability. 
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We take the inner product of 

1 d r 

2 dt 


and and sum over i, yielding 

\m^^'>\‘^dx + e [ = 0 


and hence < ||m(-^)(0)||i2(^^) < c for t < T*, where c > 0 does not depend 

on N. Next, we take the inner product of (l36|) and sum over i, and employ the 

elementary itentity (a x b) • c = (c x a) • b for a, b, c G 


(37) 


+ £ f + a 

2 dt I,. I,. 


m*'^^ X Am^^'^l'^dx 

J U) J LJ J (jJ 

= / ( — X Ho- + x x Ho-)) ■ 

J (jj 

X Ho-) • (Am*-^^ x in^^'^)dx 


-Ho 


a 

< - 
- 2 

a 

- 2 


[ |mW X AmWpdx + c [ |Hop(l + P 

J UJ J UJ 


)dx 


m 


(A) 


X Am^^^l^dx + c(l + ||m^^^||^oc(^))||H 


o-|Il°°( 0 ,T;L 2 (^)), 


where c > 0 is a generic constant independent of N. By the Gagliardo-Nirenberg inequality, 
the L°° norm of can be bounded by 

where <5 > 0 is arbitrary. The norm of can be estimated by the norms of 

Am^^) and [5l Lemma 2.1]: 

||m(^^||/^ 2 (^) < c(||Am(^^||L 2 (^) + ||m(^)||L 2 (a;)) < c(l + ||Am(^)||i 2 (^)), 

which holds for all functions with homogeneous Neumann boundary conditions. Choos¬ 
ing (5 > 0 sufficiently small, we infer from fl37jl after an integration over (0, t) with t < T* 
that 


1 

2 


|Vm*^^^(t)|^dz -I- 



jAm^^^l^dxds -|- 


a 


0 J UJ 



m 


X Am'-^^l^dxds < c{e), 


0 J UJ 


where c{e) > 0 depends on e and || Vm*^^^( 0 )|| 2 , 2 ('f^) but is independent of N. This shows 
that the solution to (l36D exists on [0,T]. Moreover, the above bound also allows us 
to perform the limit A —)■ cxd in (l36D . which gives a strong solution m G L°°(0, T; H^iuj)) fl 
L2(0, T; 7/2(0;)), dtm G L2(0, T; L‘^{u:)) to (!3^ - (|3m . 

Step 2: e-uniform estimates for m. Next, we prove some estimates for m which are 
independent of e. First, we show the bound. Let p > 1. We take the inner product of 
and |m|P“^m G -^^(0, T; L°°{u)): 

[ dxe f c}fc(|m|^~^m) • dkmdx = 0, 

duj J UJ 7_1 


d 


p+ldt 


k=l 
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where we abbreviated dk = d/dx^- Since 


^ ^m) • dkm ^ + (P “ ) dkmidkrrij > 0, 


k=l 


i,i=l k=l 

we obtain ||m(t)||ip+i((^) < ||m°||ip+i((j) for all t G [0,T] and p > 1. Exploiting the fact 
that |m°| = 1 in 01 , we may let p —)■ oo to deduce that ||m(t)|| 2 ,oo((j) < 1 for t e [0,T]. 

In order to derive some uniform gradient estimates, we take the inner product of 
and —Am and integrate in u. Similarly as in fl371) . this gives 

Id |2 , f , . ,o , <T 


(38) 


2 dt 


iVmrdx + e / |Amr(ia; + 


|m X Amrdx 


< c / |Ho-d(l -I- \iLn\'^)dx <2c \ii(^\'^dx, 


and 9/m and 


where we have used the fact that |m| < 1. Taking the inner product of 
integrating in to leads to 

[ \dtin\^dx + ~ [\Vin\‘^dx= [ 9/m • (m x (Am + 

J LU ^ J UJ J OJ 

— a 9/m ■ (m x (m x (Am + Ho-)))cix. 

J LU 

We integrate over (0, t), apply Young’s inequality, use the boundedness of m, and take into 
account estimate 



\dtm\'^dxds + - / \ViLn{t)\‘^dx — - / |Vm(0)|^(ia; 


0 J u) 


< 



\dtm\‘^ dxds + c 


< - 
- 2 


'0 J LO 

1 



|m X Am\‘^dxds + c 



'0 J (jJ 



iJiJ^dxds 


0 J u) 


\dtm\‘^ dxds + c 


0 J (jJ 



lH^rdxds. 


0 J uj 


Then, combining this estimate and the time-integrated version of (l38|) . we obtain 


(39) 


\Vm{t)\'^dx — / |Vm(0)d(i(r-I-£ 



I Am^d^ds 


+ 



|9/mddxds < c 


0 J u 



' 0 J iJ 
| 2 , 


|H„| dxds. 


0 J uj 


This gives e-uniform estimates for m in the spaces L°°{0,T; X{u)) and id^(0, T; L^(ci;)) 
and for e^/^m in L^(0, T; id^(a;)). 

Step 3: Uniqueness of solutions. Let m, m' be two strong solutions to (l3^ - (|33ll satisfying 
|m| < 1, |m'| < 1 in a; X (0,T). Set u := m — m' and recall that Ho- = (t(H -|- /5s) G 
L^(0,T; L^(a;)). Then u solves 

9/u — eAu = m X (Am -|- Ho-) — m' x (Am' -|- Ho-) 
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— a^m X (m x (Am + Hg-)) — m' x (m' x (Am' + Ho-)) 

= u X (Am + Ho-) + m' X Au — au x (m x (Am + Ho-)) 

— am' X (u X (Am + Ho-)) — am' x (m' x Au). 

Taking the inner product of this equation and u, the hrst and third terms on the right-hand 
side cancel. Then, integrating in oj leads to 

(40) f \u\‘^dx + e f iVul^dx = f u ■ (m'x Au)(ix 

dt J u) Jbj Jlo 

— a u ■ (m' X (u X (Am -|- Yl„)))dx — a u ■ (m' x (m' x Au))dx. 

J LJ J (jJ 

Integrating by parts, the hrst two integrals on the right-hand side are estimated as follows: 

/ u ■ (m' X Au)dx — a / u ■ (m' x (u x (Am -|- Ho-)))da: 

J UJ J UJ 

= — / u ■ (Vm' X '\/u)dx — / Vu ■ (m' x Vu)dx 

J LO J iO 

+ a u - (m' X (Vu x Vni))dx + a u - (Vm' x (u x Vni))dx 

J UJ J UJ 

+ a Vu • (m' X (u X Vni))dx — a / u • (m' x (u x Ho-))dx 

J UJ J UJ 

<c / |u| (I Vm|-I-I Vm'l) I Vu|dx-I-c / |up| Vm| | Vm'|da:-t-c / |up|Ho-|dx. 

J UJ J UJ J UJ 

The last integral on the right-hand side of flTOj) can be treated in a similar way: 

—a / u ■ (m' X (m' x Au))da: = a / u ■ (m' x (Vm' x Vu))dx 

J UJ J UJ 

+ a u ■ (Vm' X (m' x Vu))dx + a Vu ■ (m' x (m' x Vu))dx 

J UJ J UJ 

< c / |u| |Vm'| |Vu|(ia; — a / |m x Vu^dx. 

J UJ J UJ 

Inserting these estimates into fl40l) . we deduce that 

(41) [ |u|^dx-|-£: [ |Vu|^dx < c f |u|(|Vm|-|-|Vm'|)|Vu|dx 

dt Joj Ju) Ju) 

+ c |up|Vm| iVm'Idx-I-c / |up|Ho-|dx. 

J UJ J UJ 

The hrst integral on the right-hand side becomes 

I |u|(|Vm| -F |Vm'|)|Vu|dx < ||u||i4(^) (|| Vm||i4(^) -F || Vm'||L4(^)) || Vu||i 2 (^). 
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The Gagliardo-Nirenberg inequality and the bound |m| < 1 imply that 

||u|U4(^) < 

This shows that 

f |u|(|Vm| + |Vm'|)|Vu|da; < c||u|| (||m||^4"^) + ||m'|| 

J U) 

< <^l|u||ii(^) + c(5)||u||22(^)(||m||^2(^) + ||m'||^2(^)), 

where <5 > 0 is arbitrary. This argument is only possible in two space dimensions. Indeed, 
in three dimensions, we obtain the expression c(< 5 )||u|li 2 (^)(||in||^ 2 (^) + 11 ^ 11 ^ 2 ( 0 ;))! which 
we cannot estimate since we do not have the regularity m, m' G L'^( 0 ,T;Tr^ M). 

The second integral on the right-hand side of fl4T]) can be estimated in a similar way, 
using the continuous embedding H^{uj) 


|up|Vm| \Vm\dx < ||u||^4(^) ||Vm||i4(^)||Vm'||L4(^) 


(Also this estimate holds in two space dimensions only.) Finally, the last integral in (l4T|l 
becomes 


/ \u\‘^\H^\dx < ||u||^ 4 (^)||H^||l 2 (^) < c||u||^^i(^)||u||i 2 (^)||H<^||i 2 (^) 

J U) 

^ '^11^11^1(0;) + c(<5)||H^|lL2(a;)l|u||i2(^). 

Putting these estimates together and choosing 5 > 0 sufficiently small, we conclude from 
(ITT]) that 

\^t / +I / < c(e)(l + ^(f)) j \\i\^dx, 

where g{t) = ||H,(t)||i 2 (^) + ||m(t)||^ 2 (^) + ||m'(t)||^ 2 (^). 

As g & L^{0,T), Gronwall’s lemma and u(0) = 0 imply that u(t) = 0 in cn, f > 0, which 
finishes the proof. □ 


3.3. Uniform estimates and existence of the regularized problem. We need uni¬ 
form estimates for all fixed points of the operator F, defined in Section [2l Such an 
estimate is provided by the following lemma. Recall that X = L^(0, T; L^(G))^ and 

r = c»(|o,r];L»). 

Lemma 4 (L^ estimate). Let (p, s; m) ^ X x Y be a fixed point of F{-,a) for some 
a G [0,1]. Then there exist constants ci, C 2 (M) > 0, which are independent of e and a, 
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such that for all t G (0,T), the functional S{t), defined in ffT7|) . satisfies 



S{t) + ci / / (|Vpp + IVsp)(ix + Cl f / \dtin\'^dx < C 2 {M). 


0 Jn 



0 J (jj 


Proof. To simplify the computations, we let a = 1. The proof for general a G [0,1] is 
similar. We compute 

(42) ^ ^ + 

«/ v (jJ 

where {•, •) is the dual product between and Employing fl34l) with <7 = 1 

and p = p*, we hnd that 


(43) 


{dtp,p-pD)ds = 



DV{p - Pd) ■ (Vp - [p]M'Ei)dxds 


0 Jn 


< -- 



D\V p^dxds + 


0 Jo, 



D\V p£)ydxds 


0 Jn 


cM 



IV(p — pd)\ \Pi\dxds 


0 Jn 


< —c 



\V p^dxds + c 


0 Jn 



lEl^dzds + c(M). 


0 Jo. 


Furthermore, using fl35|) with a = 1, s* = s, and m = m*. 


(44) / {dts, s)ds = 



DlVsydxds + 


0 JQ 


J^DMmVs (|J| 


dxds 



-dxds 


iQ Jn 


< 



D\\/sfdxds + c(M) 


'0 JQ 



|Vs| |E|dxds — / ——dxds 


0 .Jn 


0 Jn 


T 


< —c 



iVs^dxds + c(M) 


0 Jn 



lErdxds. 


0 Jn 


Next, we have the energy estimate for the Maxwell equations fl28l) - fl2^ : 

D{VRl{p) ■ E - [p\M\E?)dxds - 
(|Vi?f(p)p + \dtR\{m)y)dxds + c{5) 


^ 0 J u) 




dfR^Aui) ■ Pidxds 


0 J u) 


< s 



0 j LJ 



HI dxds. 


0 J oj 
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where <5 > 0 is arbitrary. By the properties 
it follows that 
1 


and fl?rD of the regularization operators, 
/ ( 1 ^ 1 ^+ \\p(^)\\HHn)ds + 6 \\dtm{s)\\l2^^-^ds 


mi 


.^)\\L 2 (uj)ds + c(5) 



0 j UJ 


The norm of Vp can be absorbed by the corresponding term in 
sufficiently small. Estimate fl3^ can be formulated as 


^ [ \'Vm{t)\‘^dx + e f \Am\‘^dx + f \dtin\‘^dx < c /(|H 

J LlJ J UJ J UJ J UJ 

Combining the above estimates, flT2|) becomes, after time integration. 


\iiydxds. 

, choosing 5 > 0 

^ + \s\‘^)dx. 


0 Jn 


S{t) + c / / (IVpp + IVsp)(ix(is + / / \dtm\‘^dxds < c + c{M) / S{s)ds 


0 J UJ 


An application of Gronwall’s lemma ends the proof. 


□ 


Corollary 5 (Solution of the regularized problem). There exists a weak solution (p, s, E, 
H,m) to (|T9i)-(l23i) satisfying the initial and boundary conditions ([2]), ([5])-([ 6 ]), and (IloD-dllD 
as well as the regularity properties stated in TheoremUl 


Proof. Lemma 0] provides uniform estimates for the hxed-point operator F dehned in Sec¬ 
tion O Thus, the result follows from the Leray-Schauder hxed-point theorem. □ 


3.4. The limit £ ^ 0. The estimate in Lemma |4] is independent of e, which allows us to 
perform the limit e —?■ 0 . 

Lemma 6. There exists a weak solution (p, s, E, H, m) to 
(45) dtp - div(T)(Vp - [p]mE) = 0, 

Vs- [|s|]Mj^ (2) E^^ -Fym x [|s|]m|^ = 

(47) dtB - curlH = D{Vp - [p]mE), 

(48) dtH + cnil E =—dtin inQx{0,T), 

(49) dtm = m X (Am -|- H -|- (3s) — am x (m x (Am -|- H -|- (3s)) in oj x (0, oj), 

with the initial and boundary conditions (121), (ED-dSl), dH-dED, satisfying the regularity 
properties stated in TheoremU\ and the constraint |m| = 1 m ca x (0,T). 

Proof. We denote the solution to ffT^ - fl23|) with a superindex £ to indicate the dependence 
on this parameter. Lemma 0] gives the uniform estimates 

l|L2(0,r;irl(O)) + ||L2(0,T;iyl(O)) A C, 

||L2(0,T;L2(a,)) < C, 

||E‘'^^ ||loo(0,T;L2(0)) + ||H*^^^||/,oo(0,T;L2(0)) < C- 


(46) 


dfS — div 
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These estimates and show that 

ll^tP^^^llL 2 ( 0 ,T;Hh(D)') + |l^4S^^^||L2(0,r;Hl(C)') < C, 
where the constant c > 0 may depend on the truncation parameter M but not on e. There¬ 


fore, we infer from the Aubin lemma 
which are not relabeled, as e —?■ 0, 
he) 


P 


(e) 


Eh) 


^p, sh) 

—)■ s 

^p, sh) 

^ s 

^tph) - 

" dtp 

^ish) - 

^ dts 

E, Rh) - 

R 

mh) m 

mh) _ 

V * 

^ m 


dfUi 


and weak compactness that, up to subsequences 


strongly in L^(0, T; L^(r2)), 
weakly in L^(0, T; if^(r2)), 
weakly in H\){yt)'), 

weakly in L^(0, T; i7^(r2)'), 
weakly* in L°°(0, T; L^(r2)), 
strongly in C'^{[0,T]; L^{uj)), p < oo, 
weakly* in L°°{0,T; H\uj)), 
weakly in L^(0, T; L^(a;)). 


According to [221 P- 671], V-Rg(ph)) and dtRK mh)) converge weakly in to Vp and c^firi, 
respectively, taking into account (I2SD . (12711 . These convergences allow us to perform the 
limit e —0 in ([T^ - (I2211 . showing that (p, s, E, H) is a weak solution to (H7]) - (H5]) . 

It remains to pass to the limit e —)■ 0 in the regularized LLG equation fl23ll . For this, 
we observe that (I23p can be rewritten as v — amh) x v = 4, where v = mh) x (Amh) -|- 
Hh) -|- /3sh)) and 4 = c^iinh) — eAmh). The solution of this equation is v = G(amh))fg, 
where 


G(amh)) ; G(amh))f = (1 + |vp)“^ (f -|- amh) x f -|- (amh) ■ f)f]), 

is the inverse of the mapping u i-A u — amh) x u for u G Thus, fl23|) rewrites as 

G(Q!mh)) (](9jmh) — eAmh)]) = mh) x (Amh) -|- Rh) -|- /3sh)). 

Multiplying this equation with the test function 0 G C°°{u) and integrating over u, inte¬ 
gration by parts leads to 

/G(amh))0jmh)0dx-|-£ f (G(amh))Vmh) ■ V0-|-a0G'(amh))| Vmh)p])da: 

J U) J OJ 

= — I V0 • (mh) X Vmh))da: -|- [ 0mh) x (Rh) + (3s^^'^)dx. 


By the above convergence results, we can let e —)■ 0 to obtain 


G(Q!m)(9fm0(ia; = — V0 • (m x 'Vm)dx 


()m X (H -I- f3s)dx 


a.e. in (0,T), which is the weak formulation of (IT9|) . Taking the inner product of (l49ll and 
m, we conclude immediately that |m(t)| = |m(0)| = 1 in a; x (0,T). □ 
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We conclude this section by showing that (jl]) and ([7]) hold in a weak sense. Let 0 G 
C^(r2\r£)). Taking the inner product of fl47D and V0, integrating over hi, and employing 
gSD gives 

f Fi ■ V(j)dx = f curl H ■ Vcjxix + f D{Vp — [p]mE) • 'Vcjydx = f pcpdx. 
dt Jq Jq Jn dt Jq 

Consequently, 

d f d f 

(50) — / (div E — = — / E ■ 

dt Jq dt Jyj^ 

If 0 G we infer that divE(t) — p{t) is constant in time. Taking into account the 

hrst equation in (1T6|) . it follows that divE — p{t) = — C(x) holds for all t > 0 . When 
0 G C^(r2\r d), fl^ becomes {d/dt) Jp^ E ■ u(f)da = 0 and then the third equation in flTOll 
shows that E(t) ■ u' = 0 on Tat for t > 0. Finally, taking the divergence of fHHj) . it follows 
that 0tdiv(H + m) = 0 in the sense of distributions and, because of the second equation 
in flT6|) . div(H(t) + m(t)) = 0 for t > 0. 

3.5. Uniform L°° bounds for the charge and spin densities. We show that p and s 
are bounded uniformly in M which allows us to remove the truncation in fl45p - (IT7)) . 


Lemma 7. Let (p, s, E, H, m) he a weak solution to satisfying the initial and 

boundary conditions ([2]), (l5])-([n]), ffT0|) - ffT2|) . and equations dl]), ([7]), |m| = 1 in oj x (0,T). 
Then p > 0 m fix (0, T) and there exists c{T) > 0, independent of the truncation parameter 
M and the parameter (3, such that 

||p||l°°(0,T;L°°(O)) + ||s||l°°(0,T;L°°(O)) < c{T). 

Proof. Using min{0,p} as a test function in (H5|) . it follows immediately that p > 0. For 
the proof of bounds for p and s, we employ a Moser-type iteration method. For this, let 
p>2. We employ the test function [|s|]^^s/|s| G L‘^{0,T-, H^{kl)) in the weak formulation 
of fl46]) . Observing that 



= /i + h. 


The integral U is clearly nonpositive. Since 

3 / 





|s|2 J |s| 


( 52 ) 
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we can write the remaining integral Ji as 




= - f DV|s|-V([|s|]^')cia;- V [ D^^(■ Vsjdx 
+ [ D[|s|]E-V([|s|]^i)rfa:=: J 3 + /4 + / 5 . 


We show that I 3 < 0. Indeed, the dehnition of [■]m implies that V[|s|]m = X{|s|<m}V|s|, 
which yields 

h = - [ ^(P-l)[|s|]M^X{|s|<M}|V|s||^cix < 0. 

Jn 

Furthermore, since the matrix I3X3 — s (8) s/|sp is positive semidehnite, J4 < 0. For the 
hnal estimate of I 5 , we need the assumption that D is constant (also see Remark [8]). Then, 
integrating by parts and employing the hrst equation in (j4]), we obtain 


h = D{p-l) [ [|s|]^'E-V[|s|]Mda; 

Jn 

= —D [ E.V(||s|]y* = -^D f Ml,{p-C)dx 

P Jn P Jn 

>- — -D\\C\\L-o(^n) [ [|s|]^(ix > (p - 1 )F>||C||loo(q) f (j)p^M{\s\)dx. 
P Jn Jn 

Therefore, fl3T]) becomes 

^ [ (j>p,MiM)dx<ip-l)D\\C\\L--(n) [ (j)p,M{\s\)dx, 

Jn Jn 

and Gronwall’s lemma allows us to conclude that 


[\s\]l^dx<p / 0p,M(|s(-,t)|)dx 
! Jn 

<pexp ((p-l)T)||G||Loo(f^)t) f (j)p^M{\s{-, 0 )\)dx 

Jn 

= pexp ((p-l)T)||G||Loc(f^)t) f \s^\Pdx, 


In 


since |s^| < M. Taking the pth root and passing to the limit p —)■ cxo, we infer that 

s(-G)|]M||ioc(f^) < exp(T)||G||L-(Q)f)||s°||i^(f^), t > 0. 
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Now, we choose M > Mt := exp(D||C||i:,oo('f^)T)||s°||ioo(Q) and define = {x G : 

|s(x, t)| > M}. If ^Mit) has positive Lebesgue measure for some 0 < f < T, then 


M < 


M 




^ M'j' <i M, 


which is absurd. Thus, flM(t) is a set of measure zero for a.e. t G (0,T), which implies 
that |s(a:,f)| < M for a.e. x G fi, t G (0,T). Since M is arbitrary in the interval (Mr, cx)), 
we conclude that 

||s(-,t)||L°°(o) < Mt, t G (0,T). 

The proof of the boundedness of p is similar using {{[p]m — with M > K := 

max{||pD||L-(rD), l|P°ll L,00 (n)} as a test function in fIMl) (see [23]). □ 

Remark 8 (Generalizations). The boundedness result can be generalized using refined 
Moser iteration techniques. For instance, following the proof of [22], we may allow for 
nonconstant diffusion coefficients D{x) in case that the electric field E is given. It turns 
out that the L°° bounds of p and s depend on the L°°(0, T; L^(G)) norm of E. Since in our 
proof, this norm depends on the truncation parameter M, we cannot conclude the proof 
but the argument is valid if the Maxwell equations are replaced by given functions E and 
H. 

It is possible to prove that the bounds for p and s are also uniform in time. The 
idea is to exploit the gradient norm in J 3 . Using |s|P“^s as a test function in the weak 
formulation of (H 6 |) (this is possible since we already know that s is bounded locally in 
time), we find after some elementary computations that 


-/ |sF(ix + 4 
P dt Jn 


P 


-D 


P- 1 
P 


\V\s\P/^\‘^dx 


D / E-V\s\Pdx- 


T 


\s\^dx. 


Neglecting the last integral, integrating by parts in the first integral on the right-hand side, 
and employing (| 1 |), 


d_ 

dt 


\s\^dx -I- 4 


P 


-D 


P 


\V\s\P/^fdx<{p-l)D\\C\\L^^n) / \s\Mx. 

Jn 


By the Gagliardo-Nirenberg inequality, we may replace the norm of s on the right-hand 
side by its norm (by absorbing the L?' gradient norm of by the corresponding 

term on the left-hand side). This yields a sequence of recursive inequalities of the type 


dZ'p o 

-jr < ClPZp/2 + C2, 


where Zp = 


-- p ii“ilLP(0)- 

The strategy of the rest of the proof is to derive iteratively bounds for Z 2 m for all m G N, 
which are uniform in m, and to pass to the limit m ^ 00 . This can be done exactly as in 
This idea goes back to Alikakos [ 2 ]. The result is the estimate 


||s(-,t)||Lo°(o) < cmax{l, || 
where the constant c > 0 only depends on ||G||io 


^ > 0 , 


'(n)- 


□ 
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4. Free energy estimate 

We show that the relative free energy fll 8 p is nonincreasing in time under certain con¬ 
ditions. First, we comment on the spin contribution of the energy. It comes from the 
von-Neumann entropy density tr(A^logA^ — N), where “tr” is the trace of a matrix and 
N = pao -|- s • cr is the density matrix, which is a Hermitian 2x2 matrix. Here, cxo denotes 
the identity matrix and cr = (ai, < 72 , (J 3 ) is the vector of the Pauli matrices (see [30], For¬ 
mula (1)] for a dehnition). We may decompose N according to iV = p+n+ -|- p_n_, where 
p± = p±|s| are the eigenvalues of iV and n± = |(cro±(s/|s|)-cr) are the projections on the 
corresponding eigenspaces, satisfying H^ = n± and n+n„ = 0. Then, by spectral theory, 

NhgN - N = p+{\ogp+ - 1) + p_(logp_ - 1), 
which is the expression used in (IT^ . 

Proposition 9 (Monotonicity of the free energy). Let {p, s, E, H, m) be a smooth solution 
to satisfying p > |s|. Furthermore, let ||p(t)||L°o(o) < M(T), where M(T) > 0 does 

not depend on (3 but possibly onT (this is guaranteed by Lemma\^. If < 4a/(rM(T)(l-|- 
a^)) then the free energy flTS]) fulfills the ineguality 

^ + ^ [ F>(p+|Vlogp+— Ep-h p_|Vlogp_ — Ep)dx < 0, 0 < f < T. 

dt 2 Jq 

In this proposition, the diffusion constant D = D{x) is allowed to depend on x. 


Proof. We denote the von-Neumann entropy part by Espin, the electromagnetic energy by 
Eem, and the exchange energy by Egx. By computing the time derivative of Egpi^ and 
employing (IH])-(l9]) and 2p = p+ -|- p_, we hnd that 


dEgpiYi 

dt 


[ d( dtp+ log — dtp- log —^dx 
2io V PD PdJ 


IL 


D (dtp log ^ ■ dfS log —^ dx 

\ Pd |s| P-J 

1 /* ^ 

■ (Vp - pE)dx -- dJ" log — ■ 
2 Jq |s| P- 




'x -- / — log —dx 
2^0 r P- 


— Ii + h + I3 + Ia- 


The second integral becomes 

l 2 = -l f Wlog^-(V|s 
2 Jn P- 


s|E)(ia;. 
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Taking into account ([52]), we can reformulate a part of the integrand of 


3 

Zv 

i=i 


(VSj SjE) — ^ ^ i^jk 


j,k=l 


j,k=l 


SjSf^ 

jk “i pT I S. 


|S| 

3 




'^^k • E I djk 

j,k=l 


A 


2 iVSj ■ Vsfc- 


|sp / 


The matrix A = (ajk), dehned by ajk = — Sj-Sfe/Isp), is a projection and satishes 

= A. Consequently, 


3 

Zv 

i=i 


and we infer that 


2 ds ds 2 ds .n ds 


i=l 


|2 = 2 |s| 



s 


2=1 

2 


dxi dxi 


h = 


D 


s 


|s| log —dx. 
P- 


Then combining the integrals Ji and I 2 , we obtain 

1 
2 


dEsprn = f D({V\ogp+ -V log Pd) ■ (Vp+ - p+E) 


dt 


>n 


+ (V log p_ - V log pd) ■ (Vp_ - p_E)) dx 


1 


|s| log —dx. 
P- 


Next, we compute the time derivatives of Eem and 


dEp 


dt 


= / Zi)(Vp — pE) • (E — Vlogpu) — / H ■ dtmdx 

v «-/ CJ 

_ 1 

“ 2 


T> (Vp+ - p+E) ■ (E - V log Pd) 


(Vp_ — p_E) ■ (E — Vlogpo) Ida; — / H-dtmdx, 


dEg. 

dt 


= / Vm ■ VdfUidxdx = 


Am ■ diUidx. 


Adding all time derivatives, the terms involving V log pd cancel and we end up with 


( 53 ) 


dE 

dt 


1 

'2 


T)(p+|Vlogp+ — E| + p_|Vlogp_ — E| jda; 
2 ■ 


2 r 


+ 


|s| log —dx — / (H + Am) • dtuidx. 
P- Jud 
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We employ the LLG equation to reformulate the last integral: 

(54) - / (H + Am) ■ dtToa.dx 

J (jJ 

= — r(H + Am) ■ (m x (3s) + q;|(H + Am) x mp 
J iO ^ 

— q;((H + Am) x m) ■ (m x (3s)'jdx 

/ |(H + Am) X mpdx — (3 ((H + Am) x m) ■ (s — am x s)dx. 

J U) J UJ 


= —a 


At this point, we need to make some estimates. Applying Young’s inequality to the last 
integral, it follows that 


-(3 / ((H + Am) X m) ■ (s — am x s)]dx 
J U) ' 

I 


< a / |(H + Am) x m\‘^dx + 

J U) 

f \{H + Am) X m\‘^dx + ^ f (|sp + 
J LJ J LJ 


4a 

32 


s — am X srdx 


= a 


= a 


2|„|2 


a m s 


^dx 


(H + Am) X m\^dx + 13‘^—t-— 

4a 


Thus, 


becomes 


— / (H + Am) • dtmdx < (3 


2 1 T oP' 
4a 


\sPdx. 


\srdx. 


Since log((l + z)/{l — z)) > 2z for 0 < < 1, we estimate 


h = - log —dx = - ^ log 

Jn 2r P- Jn 2r 

Inserting these estimates into (|53l) . we arrive at 

dE 1 


1 + |s|/P 
1 - |s|/p 


dx < 


nzii 

Iq rp 


dx. 


11^2 - E| +p_|Vlogp_ - E| )dx 


< 




4a rp 
Since p < M{T), the result follows. 


\srdx. 


□ 
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